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Abstract. We present ROSAT HRI X-ray observations 
of all the galaxy clusters in the Le Fevre et al. arc survey 
sample in order to study the spatial distribution of the 
intra-cluster medium (ICM), and examine the expected 
number of giant luminous arcs for the sample using two 
spherically symmetric lens models constrained by our X- 
ray data. Isothermal (3 model assumes that the ICM is 
isothermal and is in the hydrostatic equilibrium. ENF98- 
NFW model assumes that the 'universal' dark matter halo 
profile proposed by Navarro, Frenk & White is a valid 
description of the underlying dark matter distribution. 
Adopting the result of N-body/gas-dynamical simulations 
by Eke, Navarro & Frenk, dark matter distribution in 
the ENF98-NFW model can be constrained by the X-ray 
surface brightness distribution of the ICM. The expected 
number of giant luminous arcs in the sample is then cal- 
culated taking into account both detection conditions in 
the arc survey and the evolution of source galaxies. We 
find that the isothermal (3 model cannot reproduce the 
observed number of giant luminous arcs even allowing un- 
certainties in the source galaxy model. The ENF98-NFW 
model displays good agreement in number of giant lumi- 
nous arcs. However, some clusters have their virial tem- 
perature 3-4 times higher than their X-ray temperature 
obtained from spectral data or from the Lx — T relation. 
Thus, we conclude that both spherical models which are 
consistent with all the available X-ray data cannot repro- 
duce the observed number of giant luminous arcs in the 
sample. To solve this discrepancy we believe that the gi- 
ant luminous arc statistics will need to take properly into 
account 'irregularity' in the mass distribution in each clus- 
ter. 

Key words: Galaxies: clusters: general - Galaxies: statis- 
tics - Cosmology: observations - ( Cosmology:) gravita- 
tional lensing - X-rays: galaxies 



1. Introduction 

A number of Giant luminous arcs (GLAs) have now been 
detected in many distant clusters of galaxies (e.g. Soucail 



et al. 1987; Lynds & Petrosian 1986; Fort & Melher 1994 



for a review) . Their spectroscopic confirmations as lensed 
normal distant galaxies (e.g. Soucail et al. 1988) lead to 



exciting new cosmological developments in this decade. 

Probing the mass distribution of distant clusters is one 
of avenues of these developments (e.g. Narayan & Bartel- 
mann |1995|; Hattori et al. |1999|; Umetsu et al. 119991). Early 



analyses pointed out discrepancy of a factor 2-3 be- 
tween cluster masses derived from strong lensing and stan- 
dard X-ray analysis of t he ICM (e.g. Loeb & Mao |l994 ; 
Miralda-Escude & Babul 1995 ). Sever al possible solutions 
are proposed: (i) Loeb & Mao 1994 proposed that non- 
thermal pressure due to equipartition magnetic field and 
turbulent may play an important role on supporting the 
ICM. (ii) A multi-phase model of the ICM in the cen- 
tral part of clusters - indicated by the cooling fiow model 
(Allen et al. |1996D 



increases by a factor 2-3 the mass 
deduced from X-ray compared to standard (single-phase) 
isothermal model (Allen 1998). (iii) Leaving from "inap- 
propriate modeling of the cluster mass distribution due to 
neglecting the contributions by substructure and member 
galaxies, which leads overestimation of the cluster mass 
from strong lensing events (e.g. Kneib et al. 1993| ; Kneib 
et al. |1995| ; Kneib et al. |1996| Hattori et al. |1998D . One can 
expect that the next generation X-ray telescopes (Chan- 
dra, XMM, Astro-E) will provide definite answer of the 
first two possible solutions. The third one seems, however, 
the closest one to the reality (Hattori et al. 



1999) 
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Instead of these detailed analyses for individual clus- 
ters, Statistics of GLAs constrains the average properties 
of cluster mass distribution. The statistics of GLAs is con- 
sist in counting the number of GLAs in a well defined clus- 
ter sample and constrains the average properties of cluster 
mass distribution. Using spherically symmetric mass dis- 
tribution models, Wu & Hammer ( 1993| ) showed that the 
predicted number of GLAs was critically dependent on 
the degree of the central mass concentration of clusters. 
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Miralda-Escude (1993a) examined arc statistics with ellip- 
tical lens models and concluded that elliptical mass dis- 
tribution did not change the above conclusion drastically. 
Using numerically simulated model for clusters, Bartel- 



2.1. Cluster sample and its X-ray observations 



and rich clusters with 



mann and his collaborators (Bartelmann & Weiss 1994; 
Bartelmann et al. 1995, henceforth BSW95) showed that 
their asymmetric cluster mass distribution (e.g. substruc- 
ture) increased the expected number of GLAs. In a further 
study, Bartelmann et al. ( |1998D examined the dependence 
of arc statistics on cosmological parameters. They con- 
cluded that the open cold dark matter model gave the 
largest number of giant arcs compared to other cosmol- 
ogy, in particular the Einstein-de Sitter universe and that 
the open cold dark matter model was likely to be the only 
model which could match current observations. 

The first well-d efined GLA survey was conducted by Le 
Fevre et al ( |l994 henceforth LF94). Hattori et al. ( |l997 , 
henceforth Paper I) proposed a method that predicted the 
number of GLAs taking into account both detection condi- 
tions in the arc survey and the evolution of source galaxies. 
They applied their method to the LF94 arc survey sam- 
ple using the Einstein X-ray data to calibrate the mass 
distribution - assuming the hydrostatic equilibrium and 
the spherical symmetry. The conclusion of Paper I is that 
the expected number of GLAs deduced from the X-ray 
analysis was inconsistent with the observed one. Hamana 



& Futamase (1997, henceforth HF97) examined the GLA 
statistics taking into account the evolution of the luminos- 
ity function found by the Canada France Redshift Survey 
(Lilly et al. 1995). They showed that the observed evo- 



lution in the galaxy luminosity function increased the ex- 
pected GLA number by a factor of 2 - 3 at most. Although 
the Lx — o relation adopted by them and by Wu & Ham- 



mer (1993) appeared to be incorrect wh en com pared with 
recen t results (e.g. Mushotzky & Scharf 1997; Markevitch 



199S ) , we can study how sensitive the GLA statistics is on 



the evolution of the galaxy luminosity function. 

This paper is the second in a series on the giant lu- 
minous arc statistics with the LF94 arc survey sample. In 
this paper we re- visit the LF94 sample with higher quality 
ROSAT HRI data because no precise measurements of the 
ICM spatial distribution in the LF94 sample clusters were 
available in Paper I. We assume that the mass distribu- 
tion in clusters are spherically symmetric. This paper is 
organized as follows. In §4, we describe the lensing prop- 
erties of mass models. We discuss the GLA statistics and 
its implication in §5, and summarize our conclusion in §6. 
Notes on individual clusters, as well as contour plot and 
radial profile of each cluster, are given in the appendix. 
Throughout this paper, we adopt a Hubble constant of 
Hq = 5O/150 km s^^ Mpc^^, the present density parame- 
ter of f2mo — 0-3, and the present cosmological parameter 
of l^Ao = 0.7. 

2. Cluster sample and its X-ray data 



4 X lO'^'^erg s^^ (Ha 



Selecting 16 distant 
ix(0.3 - 3.5keV) > 4 x iU"erg s 
50/i5o km s^i Mpc^i, O^o = 1-0, ^aq = 0.0 and 
z > 0.2) from a complete sample of distant rich clusters 
selected from the EMSS, LF94 performed a medium-deep 
V, R, and / imaging for a GLA survey and found 6 GLAs. 

X-ray observations for all the 16 clusters in the LF94 
sample were performed by ROSAT HRI. The HRI data of 
11 clusters were taken from the ROSAT Archive at Max- 
Planck- Institut fiir extraterrestrische Physik (MPE). The 
HRI data of remaining clusters were obtained by our own 
proposals. The instrument HRI had ~ 5" FWHM reso- 
lution and was sensitive for an energy range of 0.1 — 2.4 
keV (Triimper 1984 ), which makes HRI the best instru- 
ment to date to perform detailed observation of the ICM 
distribution in distant clusters. In Table |l|, we list the 
log of the ROSAT HRI observations. The column den- 
sity of the galactic hydrogen in 6th column in the Ta- 
ble is calculated using EXSAS [| command 'CALCU- 
LATE/GALACTIC_NH' which calculates the galactic hy- 
drogen column density toward the specified direction and 
is based on Dickey & Lockman (1990). Our observation 
revealed that MS 1333.3-1-1725 was not a cluster but an 
X-ray point source. 



2.2. ROSAT HRI data reduction 

We used EXSAS and XSPEC Q analysis packa ges to re- 
duce the data. The position of point sources higher than 
3cr is determined via standard source detection pipeline 
for HRI data in EXSAS. The cluster center was deter- 
mined as the brightest X-ray peak. Accuracy of ROSAT 
iffij pointing was checked with both Hambrug RASS Cat- 
alog of Optical ID (HRASSCAT) and ROSAT SIMBAD 
identifications (ROSID)^. Positions of sources higher than 
3cr were compared with positions of objects cataloged in 
HRASSCAT or ROSID. The vignetting of the ROSAT 
HRI was less than 5% within a radius of 5 arcminutes (i.e. 
600 pixels) from the detector center at all energy range 
of 0.1 - 2.4 keV for which the ROSAT HRI was sensitive. 
Thus we restricted our analysis to the inner 600 pixels of 
each image, where the background can be regarded as flat. 



2.3. Data analyses 

Photon event tables were binned into radial rings to 
make azimuthally averaged surface brightness profiles. 



^ Extended ScientiRc Analysis System provided by ROSAT 
Science Data Center at MPE 

^ provided by NASA/Goddard Space Flight Center. 

^ both are available via High Energy Astrophysics Sci- 
ence Archive Research Center (HEASARC) Online Ser- 
yice. provided by NASA/Goddard Space Flight Cent er. 
^ittp://heasarc.gsfc. nasa.gov/docs/frames/mb_w3browse. htm I 
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Table 1. The log of ROSAT HRI observations of the sample clusters. 
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Cluster name Redshift Pointing R.A. Pointing Dec. Observation date Exposure Column density 

(J2000.0) (J2000.0) [ ksec. ] [ 10^^ cm'^ ] 



MS 0015.9+1609 


0.546 


00' 


18""33=6 


+16°26 


24" 


95 Jan 


09 


- 95 


Jul 05 


76.1 


0.41 


MS 0302.7+1658 


0.426 


03' 


05'"33':6 


+17°10 


12" 


95 Aug 


19 


- 95 


Aug 27 


33.8 


1.07 


MS 0353.6-3642 


0.320 


03'' 


55"''31':2 


-36°33 


36" 


94 Jul 


19 " 


- 94 


Aug 21 


22.2 


0.12 


MS 0451.5+0250 


0.202 


04^54^09:6 


+02°55 


12" 


94 Mar 


05 


- 94 


Mar 06 


12.7 


0.78 


MS 0735.6+7421 


0.216 


07'' 


41""5a4 


+74° 15 


00" 


92 Mar 


14 


- 92 


Apr 21 


27.0 


0.41 


MS 1006.0+1202 


0.221 


10' 


08"' 45= 6 


+11°48 


00" 


96 Jun 


03 


- 96 


Dec 06 


22.4 


0.37 


MS 1008.1-1224 


0.301 


10''10™33':6 


-12°39 


36" 


94 May 


16 


- 94 


Jun 17 


68.7 


0.73 


MS 1224.7+2007 


0.327 


12'' 


27'" 14" 4 


+19°51 


00" 


94 Jun 


20 


- 96 


Jun 05 


37.1 


0.29 


MS 1333.3+1725 


0.460 


13''35™48!0 


+17°09 


36" 


97 Jul 


11 


- 97 


Jul 24 


54.4 


0.18 


MS 1358.4+6245 


0.328 


13''59™50:3 


+62°31 


12" 


91 Nov 


05 


- 93 


May 15 


29.2 


0.19 


MS 1455.0+2232 


0.259 


14''57™14=4 


+22°20 


24" 


92 Jan 


11 


- 94 


Jul 08 


14.8 


0.33 


MS 1512.4+3647 


0.372 


15''14™24^0 


+36°36 


36" 


95 Feb 


06 


" 95 


Feb 06 


35.0 


0.14 


MS 1621.5+2640 


0.426 


16'' 23"" 36^:0 


+26°33 


36" 


94 Jul 


28 


~ 94 


Jul 30 


43.8 


0.36 


MS 1910.5+6736 


0.246 


19''10™28=8 


+67°41 


24" 


96 Feb 


22 


- 96 


Mar 04 


25.1 


0.60 


MS 2053.7-0449 


0.583 


20'' 56"' 2 1^6 


-04°37 


48" 


96 May 


07 


- 96 


May 06 


4.9 


0.50 


MS 2137.3-2353 


0.313 


2l''40'''12=0 


-23°39 


36" 


92 Nov 


26 


- 94 


Apr 25 


15.5 


0.36 



The width of each ring was determined in order that the 
number of photons in each bin become greater than or 
equal to 25 to ensure that fitting could be performed, 
and that the size of each bin become greater than or equal 
to 2.5" ~ FWHM/2. A radial surface brightness profile 
was then constructed by summing up the counts in each 
bin. Note that all the contaminating point sources which 
were higher than 3cr were excluded from the photon count- 
ing. 

The radial surface brightness profile was fitted via 
Y^-minimization routine to (Cavaliere and Fusco-Femiano 



19761) 



S{0) = So 



1 



-I -3/3fit + l/2 



(1) 



where 5*0 is the central surface brightness, Of. is the angular 
core radius, and B is the background. It was physically in- 
terpreted that /3fit described the ratio of the kinetic energy 
per unit mass of the member galaxy to that in the ICM 
if cluster galaxy and the ICM distributions are isother- 
mal and galaxy velocity dispersion is isotropic. Therefore, 
the surface brightness distribution described by the Eq. 
(1) was called isothermal p model. However, it seems that 
such a situation is far from the reality (Lubin & Bahcall 
1993| ; Bahcall & Lubin |1994| ). Therefore, Eq. (1) has no 



meaning more than a conventional fitting model. We call 
the model /3 described by the Eq. (1) 'standard /3 modeV in 
this paper. The background value was first determined via 
the above fitting. Removing the background from the im- 
age we then checked that the radial source photon counts 
remained constant outside the cluster source region, if not 
we modified the background value accordingly, and we re- 
did the standard (3 model fitting. In Table ||, we list the 



standard (3 model fitting result. The bracketed numbers in 
6th column is the edge within which source photon num- 
bers are counted. 



2.4. X-ray properties of the sample clusters 

We list fiuxes, luminosities and temperatures of the sample 
clusters calculated using the best- fit values of the standard 
(3 model fitting on the left side of each column in Table |[ 
Fluxes in the 0.1 — 2.4 keV band were calculated on the 
XSPEC using Raymond-Smith thermal plasma model as- 
suming 30% of solar metallicity. The calculation of the fiux 
requires the value of the temperature. For clusters who 
have no ASCA observation, the temperature of 6 keV was 
first assumed and the temperature was calculated itera- 
tively using the 2-10 keV X-ray luminosity— temperature 
(ix — T) relation of Arnaud & Evrard (1998, henceforth 
AE98) until the temperature converges. When ASCA tem- 
perature (Tasca) was available, we used it to compute the 
X-ray flux and luminosity of the clusters. This luminos- 
ity was then used to compute the expected temperature 
from the Lx — T relation (we specify this in Table ^ by 
bracketing the temperatures with parentheses in column 



Although the Lx — T relation of AE98 was derived 
from nearby cluster sample, we used it for our distant 
cluster sample assuming no evolution of Lx — T rela- 



tion. David et al. (1993, henceforth D93) also derived 2-10 
keV Lx — T relation in (ri,„Q,17Ao)=(0: 0) cosmology us- 
ing Raymond-Smith thermal plasma model but they as- 
sumed 50 % of solar metallicity which is too high value 
for the ICM. Since AE98 showed the Lx - T relation in 
(f2,„o,f^Ao)=(l: 0) cosmology, we re-plotted the Lx — T 
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Table 2. The result of standard (3 model fitting. 



Cluster name 


5( 








I3&t 




Background 


xVd.o.f. 


count rate 








[ count s~^ 


arcsec 








[count arcsec"^ ] 




[ count s""*^ ] 


MS 


0015 


9+1609 


4.38lS;;2 


X 


10" 


6 




42.8l?-j " 
382/150 kpc 


9.88 X 10"^ t 


113.34/116 
= 0.98 


0.0344 
(< 182") 


MS 


0302 


7+1658 


6.68t^:g 


X 


10" 


6 




-I -I q-l-9.5 // 

88/150 kpc 


1.16 X lO"*^ + 


77.97/84 
= 0.93 


0.00871 
(< 221") 


MS 


0353 


6-3642 




X 


10" 


6 


o.62tE;:?» 


01 Q+18.S ;/ 
•3i.y_10.6 

208/150 kpc 


9.19to:5o X 10-'' 


112.34/113 
= 0.99 


0.0449 
(< 231") 


MS 


0451 


5+0250 


2 94+''-4i 


X 


10" 


6 


0.741°:?^ 


126.0l||:7 " 
587/150 kpc 


1.16 X 10"^ + 


100.11/110 
= 0.91 


0.158 
(< 293") 


MS 


0735 


6+7421 


5 2n+"-35 

J. ^"-0.85 


X 


10" 


5 


0.48tH? 


8.51}:^ " 

42/150"^ kpc 


1.25 X lO"** + 


159.63/117 
= 1.36 


0.0945 
(< 279") 


MS 


1006 


0+1202 


9 79+0-30 
^- '^-0.27 


X 


10" 


6 




414/i5o"^ kpc 


1.30 X lO"** + 


148.30/115 
= 1.29 


0.0437 
(< 155.5") 


MS 


1008 


1-1224 


o gi +0-76 
•-'••'^-0.61 


X 


10" 


6 


'-'•"'5-0.07 


36.4l^^«^ " 
228/150 kpc 


1.22 X lO"*^ + 


166.06/116 
= 1.43 


0.0333 
(< 183") 


MS 


1224 


7+2007 




X 


10" 


5 


'-'•*'-0.05 


6.6l*:? " 
44/150 kpc 


1.25 X lO"*^ + 


115.92/112 
= 1.04 


0.0179 

(< 170") 


MS 


1333 


3+1725 










Not 


a cluster of gal 


axies 




0.00367 


MS 


1358 


4+6245 




X 


10" 


5 


47+0-02 

'-'•^'-0.02 


rj 9 + 2.2 // 
'•^-1.6 

48/i50~^ kpc 


1.24 X 10"*^ t 


128.13/116 
= 1.10 


0.0484 
(< 212") 


MS 


1455 


0+2232 




X 


10" 


5 


'-'•"^-0.03 


69/150 ~^ kpc 


1.34 X 10"*^ + 


154.46/117 
= 1.32 


0.102 
(< 114") 


MS 


1512 


4+3647 


9 ^-,+0.58 


X 


10" 


5 




9.4t^:^ " 

68/i5o ^ kpc 


1.19 X lO"*^ + 


90.64/117 
= 0.77 


0.0178 
(< 129") 


MS 


1621 


5+2640 


o 74+0.62 
•J- '^-0.58 


X 


10" 


6 


117.1 


1246.8 " 
9739/150 kpc 


4.04 X 10"^ + 


107.36/114 
= 0.94 


0.051 
(< 166") 


MS 


1910 


5+6736 


■J-^^-0.81 


X 


10" 


6 




29.0l}« j " 
157/i5o ^ kpc 


7.20 X 10"^ + 


127.30/115 
= 1.11 


0.0167 

(< 208") 


MS 


2053 


7-0449 


< 1.1 X lO"'^ 




2/3* 




1.30 X lO"'' + 




< 0.00111 


MS 


2137 


3-2353 


±.OU_Q 20 


X 10" 


4 


63+0 04 

'-'•"'5-0.03 


8.41}:^ " 


1.22t.oSl X 10"'' 


111.01/114 


0.0728 


















54/150 ^ kpc 




= 0.97 


(< 119") 



t Fixed. See Sec. 2.4. 

J Assumed to estimate the upper limit of 5*0. 



diagram in (Omo,f^Ao) = (0-3, 0.7) cosmology. The Lx — T 
relation we used is Tx — lO'^Lx', where Tx is X-ray tem- 
perature in keV and Lx is 2-10 keV X-ray luminosity in 
erg s-\ and where (a, 6) = (-12.47 ± 0.72, 0.296 ± 0.016). 
As AE98 discussed, their Lx — T relation has the slope 
of 1/2.88 = 0.347 steeper than D93 but the Lx ~ T rela- 
tion of AE98 in (il„io,f^Ao)=(0-3, 0.7) cos mology is rather 
consistent with D93. Mushotzky & Scharf ( 1997 ) reported 
the Lx — T relation for distant clusters and showed that 
there was no evolution of Lx — T relation up to z = 0.5 
comparing with the sample of D93. In Table ||, we list 
central electron number densities (neo), central cooling 
times (icooi), ages of the universe at the cluster's redshift 
(iagc), cooling radii (^cooi), and mass-flow rates (Afcooi) for 
sample clusters. The central electron number density was 
calculated using rest frame 0.5-2.0 keV HRI luminosity 
assuming the gaunt factor of 0.9[hi'/{kBTx)]^'^'^ (Henry 



& Henriksen |1997| ), where h is the Plank constant, v is the 
frequency, ks is the Boltzmann constant, and Tx is the X- 
ray temperature, and one-ninth of the hydrogen number 
density for the helium number density is assumed. The 
cooling radius is defined as the radius where icooi = ^agc- 
The coolin g ma ss-flow rate was calculated using Eq. (2) 
in Fabian ( |1994| ). 



3. Cluster mass models and their lens properties 

In this section, we summarize lens properties of three basic 
cluster mass distribution models. 



3.1. The singular isothermal sphere (SIS) model 

Although we do not use the SIS model, the SIS model is 
employed by HF97 and then we give its lens properties 
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Table 3. The result of ENF98 /3 model fitting. 



5 



Cluster name 



So 



Background 
[ count s~^ arcsec"^ 



x'/d-o.f. 




count rate 
[ count ] 


114.33/117 = 


0.98 


0.0344 






(< 182") 


104.45/114 = 


0.92 


0.00985 






(< 221") 


116.11/115 = 


1.01 


0.0449 






zoi ) 


100.23/111 = 


0.90 


0.158 






(< 293") 


297.63/118 = 


2.52 


0.0945 






(< 279") 


149.78/116 = 


1.29 


0.0437 






(< 155.5") 


170.88/117 = 


1.46 


0.0333 






(< 183") 


130.19/113 = 


1.15 


0.0179 






(< 170") 


190.23/117 = 


1.63 


0.0484 






(< 212") 


177.84/118 = 


1.51 


0.102 






(< 114") 


97.89/118 = 


0.83 


0.0178 






(< 129") 


117.59/115 = 


1.02 


0.051 






(< 166") 


128.24/116 = 


1.11 


0.0167 






(< 208") 


132.90/115 = 


1.16 


0.0728 






(< 119") 



MS 0015.9+1609 

MS 0302.7+1658 

MS 0353.6-3642 

MS 0451.5+0250 

MS 0735.6+7421 

MS 1006.0+1202 

MS 1008.1-1224 

MS 1224.7+2007 

MS 1333.3+1725 
MS 1358.4+6245 

MS 1455.0+2232 

MS 1512.4+3647 

MS 1621.5+2640* 

MS 1910.5+6736 

MS 2053.7-0449 
MS 2137.3-2353 



0.38 

6.851^^:05 X 10--^ 
4.711°:^^ X io-« 

2.901^^9 X 10"^ 



1.38l^;l? X 10"^ 



2.83l^;2^ X 10-*^ 



3.34l°:« X 10-« 



X 10"^ 



8.30ll:^^ X 10" 



4.45;°:?^ X 10 



2.871^:^2 X 10"^ 



48.01^;° " 
429/150 kpc 

15.2tli " 
119/150 kpc 

45.61^;: " 

297/150 kpc 
135.lliy " 

629/150 kpc 

SQ.ltli " 
192/150 kpc 

70.7ltt " 
353/150"^ kpc 

52.0+41 " 
325/150 ^ kpc 

is.otii " 

119/150"^ kpc 

Not a cluster 

34.II3.I " 
226/150 kpc 

18.2tl;2 " 

102/150 kpc 

15Atli " 
111/150^^ kpc 

63.6l|g " 
497/150"^ kpc 

38.0l^:g " 
206/i5o"^ kpc 



12.7 



+1.0 // 



9.88 X lO-'^ 

LlSlHl X 10- 

9.20 X lO"'' 

1.16 X 10~^^ 

1.25 X 10"^ 

1.30 X 10"® t 

1.22 X 10"® 

1.25 X 10-® + 

of galaxies 
1.24 X 10"® 

1.34 X 10"® + 

1.19 X 10"" 
4.04 X 10"® t 

7.20 X 10"'' 

1.30 X 10"® 
1.22 X 10"® 



82/150 kpc 



t Fixed. See Sec. 2.4. 



to compare with, those of lens models we employed. The 
density profile of a cluster described by the SIS model is 



Psis(r-) 



1 



27rG r2 ' 



(2) 



where is the velocity dispersion of the cluster. The lens 
equation of the SIS model thereby becomes 



y = x- — 



(3) 



where y [x] is the angle between a source [an image] and 
the lens center in the unit of the Einstein ring radius of 



SIS model ei^^: 



qSIS 



= Att- 



D 



LS 



Dos' 



(4) 



where D^g [Dqs] is the angular diameter distance from 
the lens [the observer] to the source. 



If we assume all background sources circular, length- 
to-width ratio of a GLA equals to the ratio of the tan- 
gential stretching rate to the radial stretching rate or vise 
versa. (As one can see in the lens equation, the radial 
stretching rate of the SIS model is always unity.) The 
moduli of the eigenvalues of the Jacobian matrix of the 
lens mapping are 



At (tangential) 
Ar (radial) 





0s 




1 








61 




1 - - 

X 







dOs 



d9i 



1. 



(5) 
(6) 



and the length-to- width ratio of an arc is described as 



R: 



L 
W 



A7- 

a; 



(7) 
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Table 4. Flux, luminosity and temperature of the sample clusters. 



Cluster name Flux (0.1 — 2.4 keV) Luminosity (2 — 10 keV) Temperature^ Temperature 
[ 10-^3 g^.g g-i ^^-2 ] ^ -^QiA g-i J [ keV ] [ keV ] 



standard* ENF98** standard* ENF98** standard* ENF98** 



MS 0015.9+1609 


13.0 


40.7 


(11.1) 


o r» + 1.0 It 

8-0-1.0 


MS 0302.7+1658 


3.3 3.8 


6.7 7.7 


6.5 6.8 




MS 0353.6-3642 


16.2 


13.4 


(8.0) 


8.13+?:? 


MS 0451.5+0250 


62.1 


26.7 


9.8 




MS 0735.6+7421 


36.1 


14.8 


8.2 




MS 1006.0+1202 


16.4 


6.1 


6.3 




MS 1008.1-1224 


12.8 


11.4 


(7.6) 


7.29+?-5 


MS 1224.7+2007 


6.4 


4.3 


(5.7) 


, O+0.7 b 


MS 1333.3+1725 


0.88" 


0.12" 






MS 1358.4+6245 


17.6 


13.2 


(7.9) 


6.50+0-^ 


MS 1455.0+2232 


37.6 


17.9 


(8.7) 


5.45«:^3 


MS 1512.4+3647 


6.1 


4.3 


(5.7) 


3.57+J:? 


MS 1621.5+2640 


19.5 


38.5 


10.9 




MS 1910.5+6736 


6.2 


2.8 


5.0 




MS 2053.7-0449 


< 0.597 


< 0.659 


< 3.7 




MS 2137.3-2353 


26.4 


16.7 


(8.5) 


4.371°:^ 



* From standard P model (See Sec. 2.3) fitting results. 

From ENF98 /? model (See Sec. 3.3) fitting results, 
f From the luminosity-temperature relation in Arnaud , 

X From ASCA observations. 

i From Mushot zky fc Scharf ( |l997[ ). 
b From Henry ( |l997| ). 

^ Power low model with photon index of 1.8. 



Evrard (1995). See Sec. 2.4. 



Following Wu & Hammer (1993), a giant arc is defined 
as the image which has 

i? > eth = 10, (8) 

where eth is the threshold value of length-to-width ratio 
for identifying a giant arc (Paper I). The position of arc 
to be giant has to fulfill the following inequalities: 



< \9i\ < max{- 



eth 



7B, ——Oe) 



eth - 1 eth + 1 eth - 1 

The cross section in the source plane for forming arcs with 
R > eth hence reads 



SIS ^ I Dps^E 



3.2. The isothermal [3 model 



(10) 



Under the assumptions of spherical symmetry, isothermal- 
ity, and the hydrostatic equilibrium for the ICM distribu- 
tion, we can calculate the total gravitational mass within 
a radius r of a cluster whose ICM spatial distribution is 
described with the standard (3 model (Eq. |l]): 



M'^y{< r) = 3/3 



Gfinip/rc 1 + {r/rc)"^ ' 



(11) 



where G is the gravitational constant, fi is the mean molec- 
ular weight, TOp is the proton mass, and Vc is a physical 
core radius; Tc — DqlOc, where Dql is the angular di- 
ameter distance from the observer to the lens. We employ 
fi — 13/21. Henceforth we call this model isothermal (3 
modeF . 

The density profile of the isothermal (3 model is given 

by 

(r/r,)2 + 3 _3/3 k^T^ 1 

The lens equation of the isothermal f3 model becomes 



y ~ X — P- 



1] 



1/2 ' 



(13) 



where y [x\ is the angle between the lens center and the 
source [the image] in the unit of 9c^ and V is called the 
lens parameter: 

67r/3 /cbTx Dls 



V = 



9c firUpC^ D, 



OS 



(14) 



When P < 1 , the lens always produces a single image of a 
source. On the other hand, the lens is able to form three 
images of a source when D > 1. 
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Table 5. Central electron number density, cooling time, cooling radius and mass-flow rate of sample clusters. 



7 



Cluster name 



^cool 



^cool 





[ 10-3/150 


i/2cm-3 ] 


[ Gy 




[/i5o"'Gyr] 


[ arcsec ] 








standard* 


ENF98** 


standard* 


ENF98** 




btcLlHJ.cliLl 


IZjiN J? £70 


standard* 


ENF98** 


MS 0015.9+1609 


7.2 


6.9 


13.9 


14.5 


11.4 


U 


n 
U 








MS 0302.7+1658 


18.7 


20.2 


4.8 


4.6 


12.6 


1 A 


1 fi. 


216 


515 


MS 0353.6-3642 


7.6 


6.7 


13.3 


15.1 


13.8 


6 





4 





MS 0451.5+0250 


3.5 


3.4 


31.6 


32.6 


15.4 














MS 0735.6+7421 


40.3 


13.3 


2.5 


7.6 


15.2 


29 


35 


371 


382 


MS 1006.0+1202 


4.0 


4.3 


22.2 


20.7 


15.1 














MS 1008.1-1224 


6.8 


5.7 


14.0 


16.8 


14.1 


3 





0.4 





MS 1224.7+2007 


24.6 


17.7 


3.0 


4.1 


13.7 


18 


24 


186 


397 


MS 1333.3+1725 




















MS 1358.4+6245 


31.7 


10.4 


2.8 


8.7 


13.7 


21 


23 


348 


258 


MS 1455.0-2232 


52.5 


42.2 


1.6 


2.0 


14.6 


37 


38 


1732 


1887 


MS 1512.4+3647 


26.2 


21.6 


2.6 


3.1 


13.2 


22 


26 


475 


673 


MS 1621.5+2640 


4.2 


5.5 


27.8 


21.2 


12.6 














MS 1910.5+6736 


6.3 


5.6 


12.6 


14.1 


14.8 


12 


8 


12 


3 


MS 2053.7-0449 




















MS 2137.3-2353 


75.0 


61.3 


1.0 


1.2 


13.9 


33 


33 


2253 


2549 



★ From standard (5 model (See Sec. 2.3) fitting results. 
★* From ENF98 (3 model (See Sec. 3.3) fitting results. 



The moduli of the eigenvalues of the Jacobian matrix 
of the lens mapping are 



At (tangential) = 
Ar (radial) = 



1 - 

1 - 



V 



[a;2 + 1]V2 
V 



[x2 + l]3/2 ' 

and the length-to-width ratio of an arc is 



Ar 

Rt (tangential arc) = — 
At 



l_2?[a;2 + i]-3/2 



or 



i?r (radial arc) = — = 

Ar 



l-X)[a;2 + l]-i/2 
l_l?[a;2 + i]-i/2 



l-X>[a;2 + l]-3/2 



(15) 
(16) 

(17) 
(18) 



S.3. The universal dark rnaMer halo profile model 

NFW made series of studies on the density profiles of dark 
matter halos formed through a gravitational collapse in hi- 
erarchically clustering cold dark matter (CDM) cosmology 
using N-body simulations. They concluded that dark mat- 
ter halo reached a density profiles with a universal shape 
that did not depend on their mass ranging from dwarf 
galaxy halos to those of rich clusters, nor on the power 
spectrum of initial fluctuations, nor on the cosmological 
parameters through dissipationless hierarchical clustering. 
They found that the universal density profile can be spec- 



ified by giving two parameters: the halo mass and the halo 
characteristic (dimensionless) density 6c- 



p-NFwir) 



(19) 



Pcrit (r/r,)[l + (r/r,)]2' 

where rg is a scale radius related to the virial mass of the 
halo, and pcrit is the critical density: 



PcTitiz) = 



(20) 



&ttG 8ttG 

for Qmo + f^Ao = 1 cosmology. 

The universal dark halo profile becomes one parameter 
function in the spherical top- hat model (e.g. ENF98). The 
virial mass Mvir is defined to be the mass contained within 
the radius, Tvir, that enclose a density contrast Ac : Afvir = 
(4Ac7r/3)pcritf'vir^- This density contrast depends on the 
value of O and can be approximated by 

Ae(0, 0) = 178 0(^)°■^^ (21) 



where 



f2mo(l + -2)^ + ^^Ao' 



(22) 



for f^mQ + Q\Q = 1 cosmology. This definition shows that 
the ratio of virial radius to scale radius, which is denoted 
by the 'concentration' parameter: C = rvir/^s, uniquely 
related to 6c by 



3 ln(l + C) -C/(l + C)' 



(23) 
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Table 6. Parameters of ENF98-NFW (and MSS98-NFW*) model. 



Cluster name 5c Ts rvir Mvir Tvi 





[ Per 




[ /l50~ 


^Mpc] 


[ ^50 


^Mpc] 


[ lO^^'/iso 




[ keV ] 
ENF98 


ENF98 


MSS98* 


ENF98 


MSS98* 


ENF98 


MSS98* 


ENF98 


MSS98* 


MS 0015.9+1609 


6407 


1041 


1.14 


1.74 


5.78 


3.77 


145 


35 


34.9 


MS 0302.7+1658 


21606 


9934 


0.315 


0.400 


2.73 


2.50 


13 


8.9 


6.5 


MS 0353.6-3642 


8094 


2514 


0.787 


0.945 


4.60 


3.28 


51 


16 


15.5 


MS 0451.5+0250 


4676 


516 


1.67 


2.67 


7.92 


4.40 


214 


37 


37.6 


MS 0735.6+7421 


18019 


53857 


0.509 


0.191 


4.28 


1.61 


35 


5.5 


11.3 


MS 1006.0+1202 


5795 


804 


0.935 


1.88 


4.85 


3.88 


51 


28 


14.5 


MS 1008.1-1224 


7035 


1148 


0.861 


1.04 


4.76 


2.49 


55 


8.1 


16.0 


MS 1224.7+2007 


60568 


33733 


0.315 


0.200 


4.23 


2.12 


14 


5.1 


4.6 


MS 1333.3+1725 




















MS 1358.4+6245 


12451 


23779 


0.599 


0.218 


4.20 


2.00 


40 


3.9 


13.1 


MS 1455.0+2232 


54564 


11405 


0.270 


0.314 


3.53 


2.13 


21 


5.2 


8.2 


MS 1512.4+3647 


24594 


9266 


0.294 


0.309 


2.72 


1.89 


11 


4.1 


5.9 


MS 1621.5+2640 




















MS 1910.5+6736 


7345 


3146 


0.546 


0.714 


3.12 


2.80 


14 


10 


6.3 


MS 2053.7-0449 




















MS 2137.3-2353 


74641 


20589 


0.217 


0.245 


3.18 


2.13 


17 


5.0 


7.3 



* See Sec. 5. 



The structure of a halo mass Afvir is hence completely 
specified by a single parameter. The lens equation of the 
universal dark matter halo profile model is described as 



(Bartelmann 1996) 



y ~ X — K. 



/(a 



(24) 



where y {x) is the angle between the lens center and a 
source (an image) in the unit of the angular scale radius 
{fs/ Dol), /C is a constant coefficient: 



V ^ r DqlDlS 



OS 



(25) 



and 



Inf 



=arctanh4 / (x < 1) 

-a;2 y 1 + X \ ' 



l-ln2 {x = l) 



(26) 



Inf 



■■ arctan - 



ENF98 performed cosmological hydrodynamical and 
particle simulations to examine the evolution of X-ray 
emitting hot gas in clusters in a flat (fJmo + ^Ao = 1)j 
low-density {ilmo — 0-3) CDM cosmology. They showed 
that radial density profiles of gas in relaxed clusters are 
well described by the standard /3 model. Table 3 in ENF98 
enables us to calculate Sc and rg from the X-ray data. They 
have shown that the (3 values of most of their simulated 
clusters are around 0.79 if we exclude the cluster which is 



likely in dynamically non-equilibrium state. Parameters of 
the universal density profile is thereby specified with the 
observed core radius Tc and the normalized central gas 
densities (5gas using the relations of 



= 2.65r^ 



and 



Sr = 3.99(5, 



gas 



found in their result for /3fit = 0.79. In applying this model 
to our data, we had to remind that the best-fit /? value 
got artificially lower as the central surface brightness got 
closer to the backg round surface brightness (e.g. Bartel- 
mann & Steinmetz 1996). In general, the central surface 
brightnesses of high redshift clusters are very low and close 
to the background of ROSAT HRI. It is hence likely that 
the best- fit values of /3 shown in Table || are biased by this 
effect. To overcome this problem, we employed another /3 
model (henceforth ^ENF98 /3 modeP) fitting in which (3 
was fixed to the median value of (3 in Table 3 of ENF98; 
< Pat > median = 0.79. In this procedure, we implicitly as- 
sumed that if the radial profile was resolved up to the virial 
radius, the /3 should be ~ 0.79. We list the ENF98 /3 model 
fitting result in Table ^. We also list fluxes, luminosities 
and temperatures of the sample clusters calculated using 
the best-fit values of the ENF98 /3 model fitting on the 
right hand side of each column in Table ^. We also list 
central electron number densities, central cooling times, 
ages of the universe at the redshift of each cluster, cooling 
radii, and mass-flow rates (Mcooi) for the sample clusters 
on the right hand side in each column in Table On the 
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o 
H 




2.0 3.0 

source redshift 



Fig. 1. Total cross-section against source redshift. The 
total cross-section is calculated by summing up all the 
cross-section of each sample cluster. A cross-section for a 
sample cluster to make giant arcs of background galax- 
ies with a lens model is an area on the source plane in 
which background circular galaxies show their images with 
length-to-width ratio greater than or equal to 10 in the 
image plane. Clusters MS 1333.3-|-1725,MS 1621. 5-f 2640, 
and MS 2053.7-0449 are excluded in the calculation with 
isothermal (3 model, ENF98-NFW model, and MSS98- 
NFW model (See Sec. 5). Dashed hne: the SIS model. 
Solid line: the isothermal /3 model. Dot-dashed line: the 
ENF98-NFW model. Long-dashed hne: the MSS98-NFW 
model (See Sec. 5). 

left hand side of each column in Table ^, we list values of 
parameters of the NFW model (henceforth ENF98-NFW 
model ) derived from ROSAT HRI data using the proce- 
dure described above. The virial mass was evaluated using 
the equation (Makino et al. 1998| ) 



AnScPciitrs 



ln(l + C) 



C 



1 + C 



(27) 



For the ENF98-NFW model, one can evaluate its virial 
temperature Tvir: 



, „ 1 GMvir 
KB-'vir = -i^f^mp ■ 



(28) 



using only the information of ICM spatial distribution; 
P and 9c. Some clusters have virial temperatures which 
are much higher than temperatures obtained by ASCA or 
obtained using the Lx — T relation of AE98. 

3.4- Total cross- sections to make giant arcs 

In Fig. |l|, we show total cross-sections of LF94 sample 
to make giant arcs of background galaxies assuming that 



the background source galaxies are circular. Dashed line, 
solid line, dot-dashed line, and long-dashed line respec- 
tively represent the total-cross section to make giant arcs 
calculated with the SIS model, the isothermal (3 model, 
the ENF98-NFW model, and the MSS98-NFW model (we 
discuss this model in Sec. 5) . A cross-section of a sample 
cluster to make giant arcs of background galaxies with a 
lens model is an area on the source plane in which back- 
ground galaxies show their images with length-to-width 
ratio greater than or equal to 10 in the image plane. The 
total cross-section is calculated by summing up all the 
cross-sections of sample clusters. We calculated the to- 
tal cross-section with the SIS model using a — T relation 
used in HF97. Note that clusters MS 1333.3+1725,MS 
1621.5-1-2640, and MS 2053.7-0449 are excluded in the cal- 
culation with the isothermal f3 model, the ENF98-NFW 
model, and MSS98-NFW model (See Sec. 4.3). As shown 
in Fig. |l|, the total cross-sections calculated with the 
ENF98-NFW model and that with the SIS model are com- 
parable for all the source redshifts. Considerably larger to- 
tal cross-section of the ENF98-NFW model than that of 
the isothermal f3 model is due to their much higher Tvir. 
The main reason why HF97 could reproduce the observed 
number of GLAs in the LF94 sample is considerably high 
temperatures which they overestimated with inappropri- 
ate Lx — cr relation. 

3. 5. Properties of the giant luminous arcs 

In Table 0, observational properties of GLAs found in the 
LF94 sample are summarized as follows; major axis length: 
length-to- width ratio: l/w, distance from the cluster cen- 
ter: d, apparent magnitude: m, the projected mass within 
Einstein ring radii 6e regarding d = 9e- 



M Einstein \ 
GLA \^e) 



D 



OS 



47rG DolDls 

the projected mass deduced from the isothermal (3 model: 
37r/3 fcsTx (^ij/^c)' 



2 G^mp/r, [1-f (0b/0,)2]i/2' 



and the projected mass deduced from the ENF98-NFW 
model: 



Mf'"{eE)=Ai:p,-„,5cr,^f{xE), 



where xe = DoL^E/rs, assuming source redshift Zs — 1. 



4. The statistics of giant luminous arcs 

In this section, we examine the expected numbers of GLAs 
predicted by the isothermal (3 model and the ENF98-NFW 
model using the numerical code described in detail in Pa- 
per I. 
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4-.1. The theoretical model of galaxy evolution 



4-3. Numbers of giant luminous arcs 



As is in Paper I, prescriptions of properties of background 
galax ies are essentially the same as those in Yoshii & Taka- 
hara ( |1988D and Yoshii & Peterson ( |l99l| ). We adopted the 
galaxy type mixing ratio of (E/SO, Sab, Sbc, Scd, Sdm) = 
(0.215, 0.185, 0.160, 0.275, 0.165) given by Pence ( |1976D . 
The K and E corrections for each type were calculated us- 
ing the type-dependent, present day spectral energy distri- 
bution (SED) updated by Yoshii & Peterson ( 1991 ) and 
the type-dependent galaxy luminosity evolution models 
by Arimoto & Yoshii ( |1986| & |1987D , except for the UV 
hght of E/SO galaxies. We adopted the UV- intermediate 
NGC3379 SED (henceforth 'case I') and the UV-bright 
NGC4649 SED (henceforth 'case IF) as the SED for E/SO 
galaxies. These models are reliable, especially at low red- 
shift of z < 1 (Yoshii & Takahara |1988| , Totani et al. 
1997). Although recent observations of galaxies and their 
evolution at z > 1 (e.g. Roche et al. 1998 and references 
therein) have considerably been improving our knowledge 
of them, these models are still in good agreement with 
such recent observations; e.g. data of star formation his- 
tory taken by the Hubble Space Telescope (Madau etal. 
1996| ) or galaxy number counts (Yoshii & Peterson (1991). 



We discuss some effects of the recent new knowledge for 
the galaxy evolution (e.g. a major merger of galaxies) on 
our GLA statistics bellow. The luminosity function of all 
galaxy types was assumed to be same and was taken from 
Efstathiou, Ellis & Peterson (1988), namely the Schechter 
function with a = -1.07, 0* = 1.56 x lO-^/i^Mpc-^ and 
Mg = -19.39 -I- 51og/i {h is the Hubble constant in the 
unit of 100km s^^ Mpc^^). The galaxy formation epoch 
was assumed to he zp = 5. The absolute magnitude in 
B band for each type of galaxies was converted into V 
band magnitude using the relation My = Mb + {V — B) 
where {V - B) = -1.03 for E/SO, -0.79 for Sab, -0.64 
for Sbc, —0.56 for Scd and —0.46 for Sdm. No evolution in 
the comoving galaxy number density was assumed for the 
galaxy luminosity function. The luminosity profiles and 
the effective radii for each galaxies were modeled as the 
same way in Paper I. The intrinsic shape of the source 
galaxy image was assumed to be circular. 

4-2. The detection condition and the arc identification 
scheme 

The detection condition in the LF94 arc survey were taken 
into account as the same way in Paper I. The adopted 
GLA identification scheme was also same as described in 
section 2.4 of Paper I. The length-to-width condition and 
the apparent magnitude condition were applied for the 
images smeared by the seeing and by the limiting surface 
brightness. The threshold value for the axis ratio and the 
apparent V magnitude respectively were set to be eth = 10 
and my (arc) = 22.5. 



As noted above, the object MS 1333. 3-f 1725 is not a clus- 
ter and thus was excluded from the cluster sample. Since 
the central surface brightness of cluster MS 1621.5-1-2640 
is very close to the background. X-ray emission from only 
the central small portion was resolved by HRI (Morris et 
al. 1998). Because of this reason and poor photon statis- 



tics, the standard /3 model fitting gave unusual values as 
the best-fit result and it was impossible to estimate its 
errors. The ENF98 /3 model fitting for this cluster did not 
converge. The cluster MS 1621.5-^2640 in which a GLA 
is found hence was excluded from the arc statistics sam- 
ple in the current studies, following which the observed 
number of GLAs in the sample of 14 clusters becomes 
5. Although two giant arcs were detected by the Hubble 
Space Telescope in clusters MS 1358. 4-f 6245 (Franx et al. 



1997) and MS 1512.4-1-3647 (Seitz et al. 1997), they are not 



extremely bright and hence do not enter the GLA statis- 
tics. X-ray emission from object MS 2053.7—0449 was not 
detected in spite of pointing observation with 5 ksec ex- 
posure time. The 3(t upper limit on the X-ray flux of MS 
2053.7—0449 was set. The upper limit on the temperature 
could be obtained from this upper limit by following pro- 
cedure. First the temperature was assumed to be 6 keV 
and the upper limit of X-ray luminosity was calculated. 
Then using the Lx — T relation, the temperature value 
was reset and an upper limit of X-ray luminosity was re- 
calculated with this temperature. Iterating this procedure 
until the value of the temperature converged, the upper 
limit on temperature listed in Table ^ was obtained. The 
/3 value of MS2053. 7-0449 was assumed to be 2/3. A core 
radius of 10 arcsec for MS2053.7— 0449 was also assumed 
because smaller core radius makes a cluster more efficient 
to make GLAs. Unfortunately, we could not calculate the 
central electron number density of MS2053.7— 0449 and 
therefore MS2053. 7-0449 was excluded from the GLA 
statistics with the ENF98-NFW model (and with MSS98- 
NFW when we calculate total cross-section to make giant 
arcs). The expected number of GLAs hence becomes 4/13 
in this case. The expected numbers of GLAs in LF94 sam- 
ple, iVarc, were calculated summing up all the expected 
numbers of GLAs in each cluster, A^i.arc, which were de- 
scribed in detail in Paper I. 

We list the result in Table ||. It shows that the isother- 
mal (3 model consistent with current best X-ray data, 
cannot reproduce the observed large number of GLAs 
in the LF94's sample. On the other hand, the ENF98- 
NFW model marginally reproduces the observed number 
of GLAs in the LF94 sample. However, as one can see in 
Table ^, some clusters have extraordinary high Tvir which 
are in disagreement with the temperatures measured by 
ASCA or estimated using the L^-T relation of AE98. The 
dominant contribution to increasing the predicted number 
of GLAs is these high temperature clusters. 
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Cluster name 


I 


l/w 


d 




m 


71 f Einstein 




^NPW 




[ arcsec ] 




[ arcsec ] 


[ magnitude ] 


[IO^Mq] 


[10"Mo] 


[ 10" Mq 


MS 1006.0+1202 


5.9 


10.5 


26 


V 


= 21.92 


1.2 


0.37 


7.3 


MS 1008.1-1224 


4.5 


10 


47 


V 


= 21.53 


5.5 


1.8 


4.3 


MS 1621.5+2640 


9.8 


19.6 


16* 


V 


= 21.16 


1.0 


3.4 




MS 1910.5+6736 


6.1 


10.5 


67 


V 


= 22.29 


8.7 


1.8 


4.0 


MS 2053.7-0449 


10.5 


17.5 


15.8 


R 


= 21.91 








MS 2137.3-2353 


14.2 


17.8 


15.5 


V 


= 22.0 


0.62 


0.41 


1.1 



★ From the second brightest cluster member galaxy. For further details, see Luppino & Gioia (1992) and EUingson et al. (1997) 



Table 8. Expected number of giant luminous arcs. Case I: E/SO SED is UV-intermediate NGC3379 SED. Case II: 
E/SO SED is UV-bright NGC4649 SED. 



Model 


Case I 


Case II 


observed number 


Isothermal /3 model 


0.031 


0.055 


5* 


ENF98-NFW model 


1.7 


3.8 


4* 



* See Sec. 4.3. 



5. Discussion 

Although the ENF98-NFW model could marginally repro- 
duce the observed number of GLAs in the LF94's sample, 
some clusters have unrealistic Tvir- If ENF98-NFW model 
describes the truth and if it is this high virial temperature 
that describes the true depth of gravitational potential of 
clusters, the non-thermal pressure model (Loeb & Mao 
1994) or the cooling flow model (Allen |1998| ) would be ex- 



gives relations of 



pected. However, as shown in Table |^, cooling times of all 
the clusters with Tvir > 14 keV are longer than the ages 
of the universe at clusters' redshifts. This indicates that 
higher Tvir obtained from the ENF98 (3 model fitting were 
not due to the existence of cooling flows. The high values 
of Tvir are likely due to too large values of the scale radii 
Ts- The large values of scale radii (rg) come from large 
X-ray core radii (re) and those who have large X-ray core 
radii show signiflcant ellipticity. As NEW discussed, their 
profile came from the virialized system and those which 
have too high Tvir hence would be regarded as in relaxing 
process. 

For the comparison, we calculated a total cross-section 
to make giant arcs with another NFW type model in which 
temperatures are required to be those that are observed 
by ASCA or that are estimated using Lx — T relation of 
AE98 using the result of a theoretical work by Makino 
et al. (|199|, henceforth MSS98). MSS98 computed the X- 
ray cluster gas density distribution in hydrostatic equi- 
librium from NFW model assuming isothermality of the 
ICM. MSS98 showed that the resulting distribution was 
well approximated by the standard f3 model. Their result 



and 



= rc/0.22 - 4.55rc 
(3 = 0.96, 

27A:bTx 



and we could evaluate 6c and (henceforth we call this 
model MSS98-NFW model). We hst the parameters of 
MSS98-NFW in Table | and we plotted the total cross- 
section to make giant arcs with MSS98-NFW model with 
long-dashed line in Fig. |^. As one can see in Fig. ^ the 
total cross-section to make giant arcs is considerably small 
and almost the same as the isothermal (3 model. This 
means the NFW model which is consistent with the ICM 
spatial and spectral data of sample clusters cannot repro- 
duce the observed number of GLAs. 

This is also the same to say that what is needed is just 
to make the sample clusters' temperatures much higher to 
reproduce the observed number of GLAs. However, such 
high temperatures are no more consistent with ASCA re- 
sults or expected values from the Lx — T relation of AE98. 

Systematic errors introduced by uncertainties in the 
background galaxy model we employed are summarized 
as follows. 

As is in Paper I, t he lu minosity function was taken 
from Efstathiou et al. ( 1988 ). This luminosity function is 
in good agreement with the recent Las Campanas (Lin et 
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al. |1996| ) and the Stromro-APM (Loveday et al. |1992D red- 
shift surveys. These are called luminosity functions with 
lower normalization. On the other hand, the ESO Slice 
Project redshift survey (Zucca et al. 1997) gives higher 
normalization. Its amplitude is then higher, by a factor 
of 1.6 at M ^ M*. Luminosity functions with higher 
normalization thereby increase the number of GLAs, with 
rough estimation, by a factor of two. As noted above, we 
assumed no evolution of galaxy number density in the co- 
moving volume for the luminosity function. On the other 
hand, HF97 showed that observed evolution of the lumi- 
nosity function, whi ch ca me from Canada France Redshift 
Survey (Lilly et al. 1995 ), increases the number of GLAs 
by a factor of several. 

We investigated whether changing the type mixing 
ratio of background galaxies affected the GLA number 
by changing (E/SO, Sab, Sbc, Scd, Sdm)= (0.321, 0.281, 
0.291, 0.045, 0.062) and (0.38, 0.16, 0.25, 0.10, 0.11) and 
then little difference was found. 

As is in Paper I, the amplification factor was approxi- 
mately constant over the whole area of an image and the 
validity of this assumption has already been confirmed in 
Paper I. 

The intrinsic ellipticity of source galaxies could in- 
crease the number of GLAs by a factor of two as discussed 
in Paper L 

All these may affect on the expected number of GLAs 
by about an order of magnitude at most. Therefore 
even if uncertainties in the background galaxy model are 
taking into account, the main conclusion never change; 
the models consistent with the ICM spatial and spectral 
data of sample clusters, cannot reproduce the observed 
number of GLAs in the LF94's sample. 

As noted in Sec. 4, recent observations give us new 
insights on galaxy surface brightness and size evolution up 
to z ~ 4 (Roche et al. I998| and references therein). There 
is no size and luminosity evolution of elliptical galaxies 
at higher redshift (Roche et al. 1998| ). Spiral galaxies 
become smaller in size and brighter in surface brightness 
as redshift increases (Roche et al. 199^ ). Miralda-Escude 
( |l993b[ ) and Paper I showed that the number of GLAs 
responds sensitively on the intrinsic size of galaxies and 
the number of GLAs decreases drastically when the 
intrinsic size of galaxies becomes smaller than the seeing 



FWHM. (See Fig. 3 in Miralda-Escude 1993b or Fig. 4 in 



Paper I). We believe that this effect is stronger than that 
of surface brightness evolution because spirals seems not 
showing strong evolution in luminosity. 

If the galaxy evolution history is drawn with the 
merger model which is currently popular (e.g. Kauffmann 



bright galaxies at various redshifts. Since the current 
galaxies are to be formed by aggregation of smaller 
building blocks in the merger model, the number density 
of source galaxies at high redshift is larger than the 
current galaxy number density. These two effects may 
thus increase the number of GLAs. Although the precise 
modeling of the merger history is required to quantify the 
effect, we can regard that these effects are included as 
re-normalization in the evolution of the galaxy luminosity 
function. The result obtained by HF97 thus provided a 
rough idea how this effect changes model prediction. On 
the other hand, as we discussed above, smaller size of 
the block-building galaxies at high redshift leads drastic 
decrease of the number of GLAs. Becoming larger in size 
by merging is competed by becoming less in number of 
galaxies by merging. On the arc statistics, the effect of 
being smaller intrinsically seems stronger than both being 
more luminous intrinsically and being numerous at high 
redshift. 



1997; Bekki 1997; Bekki & Shioya 1997; Noguchi 1997) 



Asada ( 1998 ) showed that the use of the Deyer-Roeder 
distance to take into account the inhomogeneity of matter 
distribution in the universe decreased the cross-section 
for forming GLAs for all the sets of {il,A,H). Therefore 
the appropriate application of angular diameter distance 
taking account of the inhomogeneity of the universe 
further decrease the number of GLAs. 

Our calculation assumed that one GLA is generated 
from one single source galaxy. However, it may happen 
that two or more GLAs are generated from a single source 
galaxy. This means that the 'true' expected numbers of 
GLAs exists between one times the expected numbers we 
calculated and two times of them under the spherically 
symmetric mass distribution models. 

We close our discussion by this simple question: are 
all observed GLAs in LF94's sample really giant luminous 
arcs? LF94 discussed the possibility of mis-identification 
by elongated objects and claimed that 6 GLAs they 
found were really GLAs either because their widths were 
not spatially resolved, or because they presented a well- 
defined curvature. However, the giant luminous arcs in 
MS 1910.5-^6736 and MS1008. 1-1224 are located 67 and 
47 arcsec away from cluster center respectively, which are 
unusually large values that would make the cluster ex- 
tremely massive. Although, such mis-identification does 
not affect our result, it would seriously affect the GLA 
statistics which use the number of GLA on the whole sky 
extrapolating the GLA detecting rate in EMSS sample. 
Spectroscopic conformation of GLAs therefore is needed 
for all the GLAs in the sample. 



the galaxy evolution model we employed should be largely 
modified. Owing to merging-induced star formation, the 
merger model predicts the existence of temporarily very 



Kohji Molikawa et al.: The giant luminous arc Statistics. 



13 



6. Conclusions 



Appendix. 



We studied the statistics of GLAs with sphericaUy sym- 
metric lens models based on our original X-ray spatial 
data of all the 15 clusters in the LF94's arc survey sam- 
ple, obtained by ROSAT HRI. We re-investigated whether 
the observed number of GLAs in the LF94's cluster sam- 
ple could be reproduced with the cluster mass distribution 
models consistent with the X-ray data of spatial distribu- 
tion of the ICM within a frame work of spherical symmetry 
assuming isothermality and the hydrodynamical equilib- 
rium. We employed two types of cluster mass distribution 
models. One is a model comes from the conventional (3 
model used for X-ray spatial data fittings (the 'isothermal 
beta model' in this paper) . The other is the universal dark 
matter halo profile model proposed by NFW (the 'ENF98- 
NFW model' and the 'MSS98-NFW model' in this paper). 
The ENF98-NFW model is a model with the result of 
cosmological and hydrodynamical simulations to examine 
the evolution of X-ray emitting hot gas in clusters, by 
ENF98. The MSS98-NFW model is a model with the re- 
sult of the theoretical work by MSS98. Models consistent 
with current best data of spatial distribution of the ICM 
in the sample clusters taken by ROSAT HRI (isothermal 
beta model and MSS98-NFW model) gave the numbers of 
GLAs less than two orders of magnitude and this fewness 
is significant comparing to the uncertainties in the back- 
ground galaxy model we employed. On the other hand, 
ENF98-NFW model almost reproduced the observed num- 
ber of GLAs. Some clusters' virial temperatures of this 
model, however, are much higher than the temperatures 
measured by ASCA or evaluated from the Lx — T relation 
of AE98. This indicates that either: non-thermal compo- 
nents of the pressure play a significant role in supporting 
the ICM (e.g. Loeb & Mao 1994) or the mass distribution 



of sample clusters deviates significantly from the spherical 
symmetry. For example, MS1006.0-f 1202 has straight arcs 
which cannot be formed by spherically symmetric mass 
distribution. Some sample clusters show significant irreg- 
ularity in their X-ray morphologies. We believe that tak- 
ing into account the irregularity in mass distribution of 
clusters is therefore very important subject and will con- 
stitute the next step together with a better handling of 
the temperature measurement of the ICM. 
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We note down individual clusters with overlays of signifi- 
cance contour plots superposed on the Digitized Sky Sur- 
vey ^ optical images and radial profiles. A significance con- 
tour plot is made from an X-ray image smoothed with a 
Gaussian filter with a a of 15" after subtracting the back- 
ground. The lowest contour levels is 2a. In each radial pro- 
file, the solid line and the dashed line respectively repre- 
sent the result of fitting with standard (3 model and ENF98 
(3 model. The dotted line represents the point spread func- 
tion at the position of each X-ray center. 

MS 0015.9+1609 (CI 0016+16) 

One correlation was found between the brightest point 
source in the HRI field of view and objects cataloged in 
ROSID: QSO0015-hl62. As Neumann & Bohringer(|l997, 
henceforth NB97) mentioned, there are two observations 
for this cluster. Each exposure time was 70 ksec and 2 ksec. 
NB97 merged two observational data referring the position 
of the QSO. However, since l)there is only one source 
available for the pointing check of 2 dimensional position 
and 2) 70 ksec exposure time is long enough compared 
with 70-1-2 = 72 ksec exposure time, we used only the data 
of 70 ksec exposure time to avoid systematic error which 
comes from merging 2 dimensional data referring only one 
source position. Our result of standard (3 model fitting to 
HRI data is slightly different from NB97's result of the (3 
model fitting to HRI data (but consistent with it within 90 
% errors) and rather consistent with their results of that 
of PSPC data. The position of QSO0015+162 on HRI is 



^ The Digitized Sky Surveys were produced at the Space 
Telescope Science Institute under U.S. Government grant 
NAGW-2166. The images of these surveys are based on photo- 
graphic data obtained using the Oschin Schmidt Telescope on 
Palomar Mountain and the UK Schmidt Telescope. The plates 
were processed into the present compressed digital form with 
the permission of these institutions. 

The National Geographic Society — Palomar Observatory Sky 
Atlas (POSS-I) was made by the California Institute of Tech- 
nology with grants from the National Geographic Society. 
The Second Palomar Observatory Sky Survey (POSS-II) was 
made by the California Institute of Technology with funds from 
the National Science Foundation, the National Geographic So- 
ciety, the Sloan Foundation, the Samuel Oschin Foundation, 
and the Eastman Kodak Corporation. 

The Oschin Schmidt Telescope is operated by the California 
Institute of Technology and Palomar Observatory. 
The UK Schmidt Telescope was operated by the Royal Obser- 
vatory Edinburgh, with funding from the UK Science and En- 
gineering Research Council (later the UK Particle Physics and 
Astronomy Research Council), until 1988 June, and thence- 
forth by the Anglo- Australian Observatory. The blue plates of 
the southern Sky Atlas and its Equatorial Extension (together 
known as the SERC-J), as well as the Equatorial Red (ER), 
and the Second Epoch [red] Survey (SES) were all taken with 
the UK Schmidt. 
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(00'^18'"32!15, +16'i29'"25!5)(J2000). The position of X- 
ray center on HRI is (00''18™33!73, +16<i26'"07!5)(J2000) 

MS 0302.1+1658 (CI 0302+1658) 

One observation was performed for this cluster. No corre- 
lation was found between sources higher than 3 cr in the 
HRI field of view and objects cataloged in HRASSCAT 
or ROSID. The position of the X-ray center on HBI is 
(03'^05'"31!94, +17'^10'"05!0)(J2000). The position of the 
X-ray peak is consistent with that of the brightest cluster 
member galaxy within the HBI pointing accuracy. 

MS 0353.6- 3642 (S 400) 

One observation was performed for this cluster. No cor- 
relation was found between sources higher than 3 cr in 
the HRI field of view and objects cataloged in HRASS- 
CAT or ROSID. This cluster docs not show a signifi- 
cant X-ray peak, i.e. the X-ray center, we therefore de- 
fined the X-ray center as source count rate weighted mean 
position of sources around pointing center. The source 
count rates are calculated by the EXSAS command "DE- 
TECT/SOURCES" . 

MS 0451.5+0250 (Abel 520) 

One observation was performed for this cluster. No corre- 
lation was found between the sources higher than 3 cr in 
the HRI field of view and objects cataloged in HRASS- 
CAT or ROSID. This cluster docs not show a signifi- 
cant X-ray peak, i.e. the X-ray center, we therefore de- 
fine the X-ray center as source count rate weighted mean 
position of sources around pointing center. The source 
count rates were calculated by the EXSAS command "DE- 
TECT/SOURCES" . 

MS 0735.6+7421 (ZwCl 0735.7+7421) 

One observation was performed for this cluster. No corre- 
lation was found between sources higher than 3 cr in the 
HRJ field of view and objects cataloged in HRASSCAT 
or ROSID. The position of the X-ray center is consistent 
with that of the brightest cluster member galaxy within 
the HRI pointing accuracy. 

MS 1006.0+1202 (ZwCl 1006.1+1201) 

Three observations were performed for this cluster. Each 
exposure time was 3 ksec, 20 ksec, and 65 ksec. No corre- 
lation was found between sources higher than 3 cr in the 
HRI field of view and objects cataloged in HRASSCAT or 
ROSID. We used only the data of 65 ksec exposure time 
to avoid systematic error in merging the data without any 
references because 65 ksec exposure seemed long enough. 
This cluster does not show a significant X-ray peak, i.e. 



the X-ray center, we therefore defined the X-ray center 
as source count rate weighted mean position of sources 
around pointing center. The source count rates were calcu- 
lated by the EXSAS command "DETECT/SOURCES". 

MS 1008.1-1224 

One observation was performed for this cluster. No cor- 
relation was found between sources higher than 3 cr in 
the HRI field of view and objects cataloged in HRASS- 
CAT or ROSID. X-ray center (J2000.0) is (10^10™32M2, 
-12'l39™47!0). 

MS 1224.7+ 2007 

Four observations were performed for this cluster. Each 
exposure time was 14 ksec, 24 ksec, 5 ksec, and 13 ksec. 
No correlation was found between sources higher than 3 
cr in the HRI field of view and objects cataloged in both 
HRASSCAT and ROSID. One, and the only one, cata- 
loged bright point source was available for comparison of 
pointing accuracy in each image. We therefore used only 
the data of 24 ksec exposure time which was the longest 
one among them, to avoid systematic error which came 
from merging data referring only one source position. X- 
ray center (J2000.0) is (12'^27'"13!28, -|-19<^50™57!0). 

MS 1333.3+1725 

One observation was performed for this cluster. One cor- 
relation was found between sources higher than 3 cr in the 
HRI field of view and objects cataloged in HRASSCAT 
and ROSID: QS01333-I-177. 

As the contour plot and the radial profile show, this 
object is not a cluster but an X-ray point source. Therefore 
this object was excluded from the sample. 

MS 1358.4+ 6245 (ZwCl 1358.1+6245) 

Two observations ware performed for this cluster. Each 
exposure time was 14 ksec and 16 ksec. One correla- 
tion was found between sources higher than 3 cr in the 
HRI field of view and objects cataloged in HRASSCAT: 
J135903.4-F621239. However, J135903.4-F621239 locates 
near the edge of HRI field of view and was not usable 
for the reference of the pointing accuracy check. There 
was also a bright source not cataloged and thus only one 
point source near the pointing center available for the ref- 
erence of the pointing accuracy check. Since the exposure 
times of two observations were comparable to each other 
and each exposure time is not long enough, we chose the 
way to merge two files believing that pointing of each ob- 
servation is completely identical. The difference of posi- 
tion of the point source near the pointing center is 3" in 
RA and 2" in DEC and therefore we can say that these 
two positions are identical within the HRI pointing accu- 
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racy. Position of the X-ray center on HRI is (13'^59'"50!34, 
+62'^31'"04!5)(J2000). The position of the X-ray center 
is consistent with that of the brightest cluster member 
galaxy within the HRI pointing accuracy. 

MS 1455.0-2232 (IE 1455+2232) 

Three observations were performed for this cluster. Each 
exposure time was 4 ksec, 4 ksec, and 7 ksec. No corre- 
lation was found between sources higher than 3 cr in the 
HRI field of view and objects cataloged in HRASSCAT 
or ROSID. Since the exposure times of three observations 
were comparable to each other and each exposure time is 
not long enough, we chose the way to merge three files 
believing that the pointing of each observation was com- 
pletely identical. The position of the X-ray center on HRI 
is (14i^57'"14!94, -H22'i20'"35!5)(J2000). The position of 
the X-ray center is consistent with that of the brightest 
cluster member galaxy within the HRI pointing accuracy. 

MS 1512.4+3647 

One observation was performed for this cluster. Two cor- 
relations were found between sources higher than 3 a in 
the HRI field of view and objects cataloged in ROSID: 
QSO1512+370 and HD135657. The position of the X- 
ray center on HRIis (15'^14'"22!65, -h36'^36'"20!0)(J2000). 
The position of the X-ray center is consistent with that of 
the brightest cluster member galaxy within the fffil point- 
ing accuracy. 

MS 1621.5+2640 

One observation was performed for this cluster. No corre- 
lation was found between sources higher than 3 cr in the 
HRI field of view and objects cataloged in HRASSCAT or 
ROSID. The HRI image of MS 1621.54-2640 shows dou- 
ble peaks and therefore we defined the X-ray center as the 
position of the second brightest cluster member galaxy. 
(First brightest cluster member galaxy locates about 3' 
away from X-ray center (EUingson et al. 1997). ) This po- 



HRI field of view and objects cataloged in HRASSCAT 
or ROSID. The position of the X-ray center on HRI is 
(19^10'"27;94, -h67'^41™28!0)(J2000). 

MS 2053.7-0449 

One observation was performed for this cluster. No corre- 
lation was found between sources higher than 3 ct in the 
HRI field of view and objects cataloged in HRASSCAT or 
ROSID. 

In spite of 5 ksec pointing, we could have no source 
photons of MS 2053.7-0449. This indicates the low signal- 
to-noise ratio of this cluster. We gave the upper limit of the 
central surface brightness assuming /3fit = 2/3 for analyt- 
ical calculation. This cluster was excluded from our GLA 
statistics with the ENF98-NFW model and the MSS98- 
NFW model. For the GLA statistics with Isothermal f3 
model, core radius was assumed to be 10". 

MS 2137.3-2353 

Two observations were performed for this cluster. Each 
exposure time was 2 ksec and 14 ksec. No correlation was 
found between sources higher than 3 cr in the HRI field 
of view and objects cataloged in HRASSCAT or ROSID. 
Since 14 ksec exposure time is long enough comparing 
with 14 + 2 = 16 ksec exposure time, we used only the 
data of 14 ksec exposure time to avoid systematic er- 
ror which came from merging data with no references. 
The position of the X-ray center on HRI is (2lh40™15n7, 
-23'^39'"41!0)(J2000). The position of the X-ray center 
is consistent with that of the brightest cluster member 
galaxy within the HRI pointing accuracy. 
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